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Abstract 

Conserved muonic number may turn out to be a conserved muonic charge, coupled to muonic photons, 7^. Muons 
and muonic neutrinos would emit ^^'s, which might be discovered by analysing the data from the past and future 
high energy neutrino experiments (Uke CHARM II, CCFR, CHORUS, NOMAD, etc.). There are two sources 
of 7,j's in these experiments: 1) internal bremsstrahlung in pion and kaon decays into fi and Ufj,, which provide 
neutrino beams; 2) external bremsstrahlung of muons in the shielding of the neutrino beam. In both cases the "/^'s 
would pass freely through the shielding and produce narrow muonic pairs in the neutrino detectors. These pairs 
could be distinguished from the trident events (i^^ + Z Uf^ + fi'^ + ii~ + Z) hy their kinematical properties: 
small pt of the muon pair, small invariant mass, etc. All the above processes are quantitatively analysed in this 
paper. 

1 Introduction 

Speculations on the existence of massless vector particles (extra photons) coupled to baryonic number ^ 
and leptonic numbers have been under discussion for several decades. In particular, three types of such 

photons - electronic, 7^, muonic, 7^, and tauonic, 7r, - have been considered recently in ref. [||, ^ where the 
search for muonic photons at high energy neutrino experiments was discussed. In the present paper we undertake 
quantitative analysis of such a search by considering both the emission and the detection of 7;^'s. 

As it is well known, the beams of muonic neutrinos in high energy neutrino experiments are formed in a 
long decay tunnel, where pions and kaons decay : 

7r+ ^ ^+ + , TT- ^ + (1) 

K+ + K- ^ ^i- + (2) 

(K^i3 decays are less important). The decay tunnel is followed by a thick shielding of iron, earth and concrete, 
which absorbs the beam of muons and all other secondary particles (hadrons and photons). For instance, in the 
neutrino beam of the CERN SPS the 185 meters of iron were followed by 173 m of earth, and again by 20m of 
iron and 15m of concrete [^|. The beam of muonic neutrinos passes through the shielding and produces reactions 
in the neutrino detector: mainly inelastic scattering on nuclei due to charged and neutral cuiTent interactions. 

The idea of the search for muonic photons in neutrino experiments is based on the fact that the coupling 
constant of the muonic photons is small compared to a = 1/137: 

a^/a < 10-5 . (3) 

This limit follows [|[ 0| from the precision measurements of the {g — 2)^ - anomalous magnetic moment of 
the muon [^]. Because of smallness of a^, muonic photons would easily penetrate the shielding in the neutrino 
experiments and produce narrow muon pairs in the detector due to the reaction of the muon pair production by 7^ 
in the Coulomb field of a nucleus with charge Z: 

lt, + Z^ H+ +H- +Z. (4) 

On the other hand, the rates of emission of muonic photons and of production of muon pairs by them would be 
also small, which does not promise an easy hunt. 

The main background for the process (Q) are the so-called tridents: 

v^^ + Z ^ v^j_+ pi^ + p.^ + Z . (5) 
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As we will show, the dimuon invariant mass distribution for tridents is broader, and these dimuons have 
much larger pt-imbalance than in process (^, which may help to cope with the trident background. 

To get accurate numerical results for the processes discussed in this paper we used the package CompHEP 
Jiol ] to obtain analytical expressions for squared Feynman diagrams. For these purposes a new Lagrangian have 
been implemented to include new particle, muonic photon, and its interaction with muon and muonic neutrino. 
Then, for each process under the calculation, the FORTRAN code was generated for squared matrix element using 
the corresponding CompHEP option, as well as the phase space factor with proper regularization of kinematical 
singularities (see [|ll]| for the description of regularization). The nuclear form-factors were implemented in this 
code as additional factors multiplying the matrix elements of reactions and (|]). The convolution with the 
spectra of muonic photons, muons or neutrinos was realized using the so called 'structure function' option of the 
CompHEP package. The phase space integration was performed with Monte-Carlo integrator BASES, while events 
were generated by the program SPRING, both programs are combined in the package BASES/SPRING [|l2|. The 
statistical error in our MC calculations was less than 1 %. The results of CompHEP calculations are compared with 
the results of special CERN Monte-Carlo program in section ^ of this paper. 

In this paper we present the theoretical basis for the search for muonic photons in a generic high energy 
neutrino experiment. The actual analysis of experimental data is possible only for collaborations, which run the 
experiments. 



2 Emission of 7^'s in pion and kaon decays. 

The calculation of the internal bremsstrahlung in process ([l|) is straightforward, but the resulting spectrum 
has a simple analytical form only if one neglects the kinetic energy of the muon compared to (m^ — m^). Then 
for a pion decaying at rest we get: 



dB{T: fi + Vf, + 7^) 



r(7r — > iiv) 



2ap, dx 



1 2^ , 2{l-x) 
+ -x^) ■ In -) 



(6) 



where x — bj/u),nax, ^ is the energy of the muonic photon and oJmax = {m^ — '7ip/2m7r ~ rTi^r — m^. 
Expression is infrared divergent, but this divergence has no adverse physical implications, because muonic 
photons of vanishing energy cannot produce muonic pairs. Another important feature of the expression is its 
logarithmic dependence on the mass of the muonic neutrino. This logarithm is due to configurations in which the 
photon and neutrino momenta are collinear Thus, if one is lucky, one can not only discover the muonic photon, 
but to measure the mass of the muonic neutrino, even if it is extremely tiny and otherwise inaccessible. Both above 
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Table 1: Internal bremsstrahlung: differential branching dB/duj (in units of 10~^GeV~i, assuming a^/a = 10~^) 
as a function of the muonic photon energy, uj, for three different values of the muon neutrino mass, tti^^ , in the rest 
frame of decaying meson. The maximum photon energy is comax = ("^^/k ~ "^m^/^^it/K- 
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Figure 2: Comparison of the spectrum of muonic pho- 
Figure 1: Spectra of muonic photons from internal tons from internal bremsstrahlung in the CERN 450 
bremssti-ahlung in the CERN 450 GeV neutrino WBB GeV WBB neutrino beam for two masses of muon neu- 
(Wide Band Beam). The muon neutrino mass is as- trino (dotted Une m„ ~ 10~^" eV and dashed line 
sumedtobe rrii, = 10^^" eV. — 100 KeV ). The spectrum of muonic photons 

from external bremsstrahlung is shown by a solid line. 



features (infrared and coUinear) are characteristic for the spectrum of internal bremsstrahlung in a general case, 
when the kinetic energy of muon is taken into account. To obtain accurate numerical results for the differential 
branching, in particular in the case of K decay, we evaluated analytically the squared matrix element with the help 
of the package CompHEP Then, the integration over squared invariant mass of neutrino and 7^ was done 
analytically using the computer algebra system REDUCE Jl^]. The obtained expression was evaluated numerically 
with high accuracy by the corresponding REDUCE option. The results of these accurate numerical calculations, 
differ from (||) by a few percent in the case of pion. For the pion and kaon decays they are presented in table [l| for 
a set of neutrino masses, assuming — lO^'^a. 

A muonic photon with energy uj in the rest frame of a meson (pion or kaon) and angle 9 (in the same frame) 
with respect to the direction of the meson beam has in the laboratory frame energy E^^^ , where 

E^^ = LJ7(1 + i;cos6') . (7) 

Here v is the velocity of the meson and 7 = l/Vl — t'^. The spectrum of muonic photons is derived by integrating 
over decs 9 and convoluting with the spectra of decayed pions and kaons. Photons are emitted within a cone with 
an angle 9' ~ I/7. Thus, one should take into account the angular acceptance of the detector determined by the 
geometry of experiment (the length of the decay tunnel and the position and cross dimension of the detector). In 
this paper we use as an example the CHARM II neutrino detector at the 450 GeV CERN SPS neutrino beam. 
This detector is located at 882 m from the tt/K production target with the length of the decay tunnel of 414 m. 
The formation of the 7^ spectrum is discussed in [ [l4| ] in details in the framework of a simplified approach to the 
simulation of the detector geometry and meson beam. As an example, we used the meson spectra from ref. with 
the effective K/7r-ratio equal 0.24. ^ 

The resulting spectra of muonic photons, produced by the internal bremsstrahlung, are presented in fig. |l]. 

An important feature of fig. [l| is that hard muonic photons are mostly produced in the kaon decays. Of course, 
it depends on geometrical parameters of the experiment (see [[l4|]). However this feature could be taken into ac- 
count in future experiments to increase possible yield of muonic photons produced by the internal bremsstrahlung. 

In order to see the dependence of rates and spectra on the mass of the muon neutrino we present in fig. ^the 
spectra of muonic photons for two values of the neutrino mass: m^^ ~ 100 KeV - close to the experimental limit 

We made the calculations only for the neutrino component of the neutrino beam (focused positive mesons). Similar calculations should 
be done for the antineutrino admixture due to wrong sign mesons and for antineutrino beam (negative mesons focused), where we expect 
similar results. Also we note that in this semi-analytical approach we neglect angular and radial distributions of mesons and assume the 
uniform distribution of the meson decay point along the decay tube. The exact account for these effects and their correlations is possible by 
pure Monte-Carlo simulations and is presented in section rj. 
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and another assumed to be very small, rrijj^ = 10^^" eV. The yield of muonic photons is about seven times 
higher for the second choice of . This allows one to hope that the mass of muonic neutrino could be measured 
in the experiments discussed, if muonic photon is discovered. 



Emission of 7^'s by muons in the shielding. 

The effective material for the muon bremsstrahlung in the shielding is iron; 



Fe 



M + 7m + 



Fe. 



(8) 



The cross-section of the bremsstrahlung for the electron has been calculated by Bethe and Heitler [T^]; 
its derivation may be found in many text-books, e.g. [ p^ [T^. In the case of muon there is an obvious substitution 
me m^. Moreover, the nuclear form-factor must be accounted for. This has been done by a group at the 



Moscow Engineering Physics Institute (MEPhI) 



■[E2|, see also 



In order to apply their expressions to 



the process we have to multiply them by a^/a. In the ultra-relativistic limit, when Ei ^ and E2 3> 
where Ei {E2) is the energy of the initial (final) muon, we get: 
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(9) 



Here E-y^ is muonic photon energy, Z is the charge of a nucleus of iron [Z = 26). The term 1.5 appears due to the 



nuclear form-factor, which was assumed in the form^: 

F{q^) = exp(-q2a2/6), 



(0.58 + 0.82Ai/3)fm. 

2|, 



Another parametrization of the nuclear form-factor correction is suggested in ref 
of 1.5 in the second bracket of eq. (m, one has to substitute 



(10) 

according to which, instead 



Dn 



where £>„ 



2)/m^ 

56 for iron), 5 



(11) 



mlEy/2EiiEi-Ey) 



1. 54 A"- A - mass number of the nucleus (A ^, ^ ^ ^^^^ 

minimal momentum transfer, and e = 2.718. The final numerical results of eqs. (|) and ^ are close to each 
other. In this paper we do not take into account the atomic screening and inelastic form-factors, both nuclear and 
atomic (see ref. [^]). 

The total cross-section of the process (||) convoluted with the spectrum of muons emitted in pion and kaon 
decays equals to 1760 pb for = 10~^a. This result was obtained with the help of programs CompHEP 
and BASES. To calculate the spectrum of muonic photons emitted by muons in the iron shielding we used the 
formula (^ convoluted with the spectrum of muons. The spectra of muonic photons, produced by the external 
bremsstrahlung in the iron shielding are presented in fig. ||. The internal and external spectra are compared in 
fig. H First of all, one can see that the yield of muonic photons produced by the internal bremsstrahlung is higher 
for mi/, = lO^^'^eV than the yield by the external bremsstrahlung, supporting prospects to measure muonic 
neutrino mass if 7^ is discovered. The internal bremsstrahlung spectra are harder than the external ones. Thus in- 
creasing the energy of mesons (especially of kaons) should help in measuring neutrino mass. Moreover, the ratio of 
pion and kaon components is drastically different for muons (and hence for their bremsstrahlung in the shielding) 
from that for the mesons themselves. Indeed, while initial meson spectrum in the CERN beam has 80.4% decayed 
pions, the muon flux at the end of the decay tunnel has 96.5% muons emitted in the pion decays. The relative 
enhancement of pion component is caused by a strong suppression of the kaon component, because the muons and 
muonic gammas from K-decays have wider angular distribution and are partly missing the acceptance of the neu- 
trino detector (see [|l4|] for the corresponding formulas). However the contribution of energetic muons from kaon 
decays is enhanced by rather high energy cut, 7 GeV, on muonic photons and the increase of the bremsstrahlung 
cross-section with energy. Finally, 82.9% of muonic photons and 80.3% of muon pairs originate from pions (this 
effect is discussed in more details in ref. ||l4|]). 

In the external bremsstrahlung calculations we have neglected angular distribution of muonic photons in 
the reaction (^. Indeed, muons, moving at angles larger than the angular size of the detector could, nevertheless, 
produce muonic photons crossing the detector Moreover, muons, moving in the directions crossing the detector, 
could produce muonic photons going outside it. We assume that these two effects more or less compensate each 
other. The angular distribution of muonic photons produced in the external bremsstrahlung has a forward peak 
with characteristic polar angle 6 ^ nif^/Ei, where Ei is the muon average energy. For = 31 GeV , more than 
80% of photons lie inside this 6'-cone (about 0.2°). Thus, we assume that the effect under discussion should give a 
correction factor not larger than 2 in our final results. 

An extensive compilation of nuclear elastic form-factors can be found in ref. fcffl. 
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Figure 3: Spectra of muonic photons from external 
bremssti-ahlung in the CERN 450 GeV WBB neu- 
trino beam. The contribution from 7r-decay is shown 
as dashed line, from K-decays by dotted line. 



Figure 4: Invariant mass distributions (in units of pb/20 
MeV) of muon pair from the conversion in the neu- 
trino detector of muon gammas with fixed energies of 
20 GeV (solid), 50 GeV (dashed) and 100 GeV (dotted 
line). 



4 Photoproduction of muon pairs by 

The photoproduction of muons by ordinary photons has been analysed in a number of papers. The corre- 
sponding cross-section is connected with the bremsstrahlung cross-section by crossing (see refs. [|l7[ Q). One can 
use the differential cross-section suggested in ref. [|7|], substituting instead of a, but keeping the factor Za 



Z) = 



2Za\j 



In ; 



2 m 



1 - 



A 



3 me ^2/3 ) 



3 

In ( 1 



dx+ X 

me 



(12) 



Here a — 1/137 , Z is atomic number of the nucleus {Z = 10 for the 'average nucleus' in the CHARM II 
detector), 6 = m^/2£'^^a;+a;_ - minimal momentum transfer, E^^ - the energy of the muonic photon, x± 



E±/E-y^, where E± - energy of /i 
logarithm. 



^0 



190 - a constant, which determines the value of the radiation 



With£;^ r 
factor simplifies: 



20 GeV and 0.1 < x± < 0.9, the atomic screening may be neglected, so that the logarithmic 
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Another, more recent parametrisation of the same cross section^ is 
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where B = 183, £»„ = 1.54 • A^-'^'^ and A is a number of nucleons in a nucleus (A = 20 for the 'average nucleus' 
of the CHARM II detector). Note that for not too small values of the S parameter (S/rrie ^ Z^l'^ j B-,/e) the square 
brackets in eq. ( [l4| ) reminds the last factor of eq. (^ with substitution of the term (—1.5) by eq. (jll]): 



{} 



In 



1 Du 

2 ^ l + 5{D^^-2)/m^ 



(15) 



^) S.R.Kelner, R.P.Kokoulin, private communication 
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Figure 5: Transverse momentum distributions (in units 
of pb/3 MeV) of muon pair from the conversion in the 
neutrino detector of muon gammas with fixed energies 
of 20 GeV (soHd), 50 GeV (dashed) and 100 GeV (dot- 
ted Une). 



Figure 6: Azimuthal angle distributions (in units of 
pb/1.8°) between two muons from the conversion in 
the neutrino detector of muon gammas with fixed ener- 
gies of 20 GeV (soHd), 50 GeV (dashed) and 100 GeV 
(dotted line). 



Note also that if we neglect both atomic screening and nuclear form-factor and integrate over dx+, we reproduce 
an expression similar to the well known result [ref. |l3], eq. 94.6]: 

7 /2ZaV / 2E^^ 109 \ 
'' = ^^9[-^) ['^^^^)- '''' 

Similarly to the case of the muon bremsstrahlung in the shielding the nuclear elastic form-factor can be accounted 
for effectively by introducing an additional term in brackets, —1.2. This value can be obtained from the parametri- 
sation (|l^ for average nuclei with Z = 10, 'neon', in the CHARM 11 detector Note that this term equals to —1.5 
for Fe (Z = 26), see (|). 

To get accurate numerical results for cross-section and to calculate distributions for dimuon photoproduc- 
tion we used the programs CompHEP and BASES/SPRING . In particular, the total cross-section of dimuon pho- 
toproduction by an "average" muonic photon, for — lO^^a, is 41.5 pb in the case of convolution with the 
normalized spectrum of muonic photons produced by the external bremsstrahlung, and 35.6 pb in the case of the 
internal bremsstrahlung (with m,y^ = 10^^" eV). Note that in the case of internal bremsstrahlung the shapes of all 
distributions depend only slightly on the neutrino mass. Of course this is caused by the similarity of muonic photon 
spectra for different neutrino masses (see fig. This similarity can be easily understood as the main contribution 
comes from the logarithmic terms like that in eq. (^. 

Returning to the formula ( p^ we note that its modified version, with additional term —1.2 effectively ac- 
counting for elastic nuclear form-factor, differs from the accurate calculation by CompHEP and BASES by less 
than 1% for E^^ > 50 GeV, and this difference does not exceed 10% for lower energies within the energy range 
Ej^ > 7 GeV discussed here. 

Muons should have high enough energy to be detected; in our calculations we took the detector energy 
threshold for muons as 4 GeV. Keeping in mind this energy threshold we have applied energy cut 7 GeV on the 
energies of the original muons and neutrinos and on the energy of muonic photons. However for the final muons, 
produced in the detector by muonic photons, we have applied cut E"^ > 4 GeV. All values of the photoproduction 
cross-section given above were obtained with these cuts. 

Taking into account the yields of muonic photons, normalized to one muon neutrino from two-body tt/K 
mesons decays, calculated as discussed above (see [p^), 

3.1 • 10^^ ext. brems., 

1.5 • lO"'^ int. brems., m^^ = 100 KeV, 

10.4 • 10"'^ int. brems., mj^ = 10"i° eV, 

one can derive the the specific yields of muon pairs per unit neutrino flux: 
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Figure 7: Energy spectra of muons from the conversion 
of muonic photons and from tridents in the neutrino de- 
tector at CERN neutrino beam. Solid Une corresponds 
to external brems, dashed line to internal brems and 
dotted line to tridents. 



Figure 8: Energy spectra of muon pairs from the con- 
version of muonic photons and from tridents in the neu- 
trino detector at CERN neutrino beam. Solid line corre- 
sponds to external brems, dashed line to internal brems 
and dotted line to tridents. 



2.1 • 10^^*^ ext. brems. events cm^ /neutrino, 

0.9 • 10"^*^ int. brems. events cm^ /neutrino , ^ — 100 KeV, 

5.9 • 10"^*^ int. brems. events cm^ /neutrino, m^^ = 10"^" eV. 

These specific yields of muon pairs are not exactly proportional to the specific yields of muonic photons because 
of the difference of the mean energies and effective cross-sections. To obtain these numerical results we took 
the number of average nuclei (with Z = 10, 'neon') in the CHARM II detector equal to 0.16 • 10'^^ [||, and 
= 10~^q;. Neutrino flux in CHARM II experiment is estimated [ ^8[ ] as 8.6 • 10^^ v/cnP integrated over the 
whole energy spectrum. So we would expect 177 dimuon pairs in CHARM II exposure from muonic photons 
produced by external bremsstrahlung. Then, in the case of internal bremsstrahlung we should observe about 74 
events, if the mass of muonic neutrino is close to the experimental limit, 100 KeV, and more than 511 events, if 
the mass of muonic neutrino is very small, of order 10^^" eV. Invariant mass and transverse momentum of muon 
pairs, azimuthal angle between two muons distributions are shown on figs. ^, ||, ||at some fixed energies of muonic 
photons: 20 GeV (soHd lines), 50 GeV (dashed lines) and 100 GeV (dotted fines). 

The shapes of the transverse momentum and azimuthal angle distributions scales with energy, while the 
invariant mass of two muons slightly increases at high energy. 

The energy spectra of muons and muon pairs from tridents and from conversion of muonic photons at CERN 
neutrino WBB are presented in figs. ^ and ^. 

In these and following figures we represent external brems by solid lines, internal brems with = 10~^° 
eV by dashed lines and tridents by dotted lines. All distributions are normalized to one trident event. One can note 
that these spectra are harder in the case of the internal bremsstrahlung than in the external one. 

One can see on fig. ^ that muons from 7^ conversions are produced at small angles: more than 80% at angles 
less than 1° with respect to the direction of generic muonic photon. We have found also that in both cases (external 
and internal ones) muons are produced with rather small pt: 90% of muons have transversal momentum less than 
200 MeV (fig. |l^). 

We also calculated distributions in different kinematical variables to find most effective cuts to localize the 
effect of muonic photons. In figs. [l^, |l4| we present the distributions in the invariant mass of dimuons, in 

the transversal momentum of a muon pair (the so-called -imbalance), in the azimuth angle between muons (it is 
an angle between the transversal components of and /Lt~ momenta), and in the angle between total momenta of 
fi^ and /i^. 

5 The trident background. 

As it has been mentioned in the introduction, the main background for the search of muonic pairs produced 
by muonic photons is presented by tridents, process (^. The cross-section of trident production in a Coulomb field 
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Figure 9: Angular distribution of muons from the con- 
version of muonic photons and from tridents in the neu- 
trino detector at CERN neutrino beam. Solid line corre- 
sponds to external brems, dashed line to internal brems 
and dotted line to tridents. 



Figure 10: Transverse momentum of muons from the 
conversion of muonic photons and from tridents in the 
neutrino detector at CERN neutrino beam. Solid line 
corresponds to external brems, dashed line to internal 
brems and dotted Une to tridents. 




Figure 1 1 : Invariant mass of muons from the conver- 
sion of muonic photons and from tridents in the neu- 
trino detector at CERN neutrino beam. Solid line corre- 
sponds to external brems, dashed Une to internal brems 
and dotted line to tridents. 



Figure 12: Transverse momentum of muon pair from 
the conversion of muonic photons and from tridents in 
the neutrino detector at CERN neutrino beam. Solid 
Une corresponds to external brems, dashed Une to in- 
ternal brems and dotted Une to tridents. 
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Figure 13: Azimuthal angle between muons from the 
conversion of muonic photons and from tridents in the 
neutrino detector at CERN neutrino beam. Solid line 
corresponds to external brems, dashed line to internal 
brems and dotted line to tridents. 



Figure 14: Angle between muons from the conversion 
of muonic photons and from tridents in the neutrino de- 
tector at CERN neutrino beam. Solid line corresponds 
to external brems, dashed line to internal brems and 
dotted line to tridents. 



of a nucleus with charge Z, can be easily estimated using Weizsaecker- Williams approximation (see e.g. [p9[]): 

/ I v I + I - I '7\ 2(Zq;) G , 2 , 2 n 1 ^max 

Cr(j/^ + Z ^ Z/^ + + ^ +^) = -^j, 9r) ■ Smax-y^^^ ■ (17) 

Here G — 1.166 • 10^^ GeV^^ is the four-fermion coupling constant, gh — \ + sve? 6w, Qr — sin^ 6w, 
sin^ 9w — 0.23, Smax = '^EiQmax, Qmax = l/o, wherc a is the radius of the nucleus, see eq. (p^; Ei is 
the energy of the initial neutrino. This approximate formula gives the value smaller by ^ 10% then the accurate 
calculation using CompHEP and BASES within the discussed energy range. The cross-section convoluted with the 
normalized neutrino spectrum gives the value a = 0.045 • 10"'' pb. This numerical result has been obtained with 



the nuclear form-factor (10) and with energy cut > A GeV applied on muons. The number of trident events is 
estimated by using the neutrino flux and target mass given in Section ^to be 50 events. 

Energy spectra of muons, and muon pairs for tridents are presented by dotted lines in figs. ^ and ^ In 
figs, n , 12[ |l^ |l3| we present distributions (also by dotted lines) in different variables for tridents — in invariant 
mass of the muon pair, in transversal momentum of the muon pair (pt-disbalance), in the angle between muon and 
antimuon, and in azimuth angle between muon and antimuon. One can see that trident events are distributed in 
these variables broader than muon pairs in dimuon photoproduction. The characteristic values of kinematical cuts 
on these variables are collected in table ||, together with the percentage of the background trident events which 
survive these cuts. It is obvious that the most promising is the cut in the transversal momentum of muon pair, 
P4-disbalance. Only 2.12% of trident events have p^^ < 70 MeV. By applying this cut one can hope to get the 
suppression of trident background by a factor 50 while practically all muon pairs from photoproduction will be 
kept. 

The cuts on invariant mass of muon pair and on azimuth angle between two muons could be also effective, 
each could suppress trident background by a factor 3-4. In addition, we found that distribution in transversal 
momentum of muons is much broader for tridents than for dimuon photoproduction, here the suppression factor 4 
could be also obtained for tridents. Of course special analysis is necessary to find those combinations of these cuts 
which will give maximal suppression factor for the trident background - the correlations are important. 



6 Monte-Carlo simulation of the production and detection of the muonic photons. 

In order to take into account the details of the geometry and kinematics of the neutrino beam line, a special 
CERN Monte-Carlo program was used, which directly simulates the decay of tt/K with muonic photon and the 
bremsstrahlung of muonic photons in the iron shielding of the neutrino beam. This program uses the Monte-Carlo 
vr/K-decays produced by GBEAM program [|]| with realistic geometry, magnetic focusing in the horn and reflector 
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Observable 


int.brems. 


ext.brems. 


tridents 


M^^(MeV) < 


500 


500 


39% 


p^''(MeV) < 


72 


73 


2.12% 


p^(MeV) < 


200 


200 


24.9% 


fnnideg) > 


142 


142 


25.6% 


z9^^(deg) < 


2.8 


2.7 


64.2% 


^,,(GeV) < 


28 


28 


65.8% 


E,,^{GeW) < 


56 


52 


69% 



Table 2: The kinematical cuts, each of them preserves separately 90% of muon pairs produced by muonic photons 
in the detector (columns 2 and 3), and the percentage of the trident background events within the cuts given for 
external bremstrahlung case (column 4). Here A/^^ is the invariant mass of a muon pair; and p^^ are transversal 
momenta of a muon and a muon pair correspondingly, -d^^ is the angle between and while (^^^ is the angle 
between their transversal momenta; i?^ - energy of a muon (/i+ or and is energy of a muon pair For 
internal bremsstrahlung the numerical results only slightly depend on neutrino mass within the discussed range, at 
least up to 10"^" eV. 

and tt/K lifetimes. So at the input every event has the coordinates of the vr/K decay point and the vector of tt/K 
momentum. The difference with the semi-analytical approach is that the radial and longitidinal spatial distributions 
of tt/K in the decay tunnel, the angular distribution of tt/K mesons and their correlations are automatically taken 
into account. Both positive and negative tt/K mesons are used for muonic photon flux calculation in the neutrino 
Wide Band Beam (WBB), where positive mesons are focussed, while negative mesons are defocussed. 



Observable 


Semi-analytical 


Monte-Carlo 




16.8 GeV 


17.9 GeV 


^ext 


18.5 GeV 


22.8 GeV 


10.4- IQ-^ 


12.5- 10"^ 


F^^'/F, 


3.1 • IQ-'^ 


2.3 • 10"'^ 


—int 


35.6 pb 


39.9 pb 


—ext 


41.5 pb 


53.8 pb 


jyj-tnt 


511 events 


688 events 




177 events 


172 events 


Ajtridents 
total 


140 events 


138 events 


Ajtridents 


49 events 


47 events 



Table 3: Results of semi-analytical and Monte-Carlo calculations of the following quantities: the average energies 
of muonic photons (i?*"*, E";^^*) , the ratios of muonic photon flux to neutrino flux for external and 
internal F!^^^ / F^ bremsstrahlung, the average muonic gamma interaction cross-sections for gammas from internal 
(T™* and external cr^^* bremsstrahlung, the numbers of events of muonic gamma conversions from internal TV'"* 
and external A^^^* bremstrahlung and the numbers of trident events without (A^jotaf"*") and with (^*QgY**j) the 
cut on muon momentum of 4 GeV. The fluxes, cross-sections and mean energies were calculated with the cut 
E^>1 GeV and > 4 GeV. 

The internal bremsstrahlung of muonic gammas in tt/K decays is generated according to the formulas used 
in the section |[ The flux of the muonic photons is calculated in the fiducial volume of CHARM-11 detector with 
square cross-section of 3.2x3.2 m^ located at 882 m from the proton target. 

The production of the muonic photons by external bremsstrahlung from muons, stopping in the iron shielding 
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was calculated according to the formula (^. For this study the shielding was appriximated by 80 blocks of iron, 
each 3m in length. The simplified tracking take into account the energy loss of 4.5 GeV per block and multiple 
scattering of muons in the iron. The muons were tracked down to 10 GeV. Muonic photons were generated with 
the angular distribution (f.w.h.m. 0.12°) around the muon direction at each tracking step. The effect of the multiple 
scattering reduce the flux of muonic photons in the volume of the neutrino detector, because of the increase of the 
angular spread. 

The conversion of muonic photons into muon pairs in the target of the neutrino detector was simulated 
according to formula (|l^ for the average CHARM II nuclei (A=20, Z=10). The energy cut of 4 GeV was applied 
on both muons. 

The cross-section and rates of the trident background were calculated by Monte-Carlo program |^] with 
neutrino spectra predicted by GBEAM ^ program. Only coherent part of the cross-section was estimated for this 
study. 

The results of these calculations in the form of the mean energies, cross-sections and rates are summarised in 
the table |]for the neutrino mass m^^ = lO^^*' eV and a^/a = 10^^. The numerical results of this Monte-Carlo 
simulation are close to the semi-analytical calculations. The differences are due to different approximations and 
assumptions on the geometry, the spectra of mesons, the effects of multiple scattering and the angular distributions. 

7 Conclusions. 

Let us summarize the main results of this paper If ^ lO^^a, as allowed by {g — 2)^ data, then in the 
high-energy neutrino experiment of the type of CHARM II one may observe ^ 180 narrow muon pairs produced 
by muonic photons in the case of external bremsstrahlung. In the case of internal bremsstrahlung there may observe 
about 70 events if the mass of muonic neutrino is 100 KeV, and ~ 600 events if the mass of muonic neutrino is very 
small, of order 10^^° eV. This should be compared with approximately 50 muon pairs from trident events. Thus, 
muonic photons could be discovered in such experiments if ~ lO^^a. Moreover, one may hope to measure the 
mass of muonic neutrino if this mass is small enough and muonic photon is discovered. 

By selecting pairs with pt < 70 Me V one may reduce the number of background trident pairs by a factor of 
50, to a few events. One may hope to further suppress the trident background by selecting muon pairs with invariant 
mass less than 500 MeV, or with azimuth angle between muons larger than 140°, or with transversal momentum 
of muon and antimuon smaller than 200 MeV. Thus, if no events satisfying the cuts are observed, one can reach 
the sensitivity of a^/a < 10~^. 
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